Large vortex state in ferromagnetic disks. 
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Magnetic vortices in soft ferromagnetic nano-disks have been extensively studied for at least sev- 
eral decades both for their fundamental (as a "live" macroscopic realization of a field theory model 
of an elementary particle) as well as applied value for high-speed high-density power-independent 
information storage. Here it is shown that there is another vortex state in nano-scale ferromagnetic 
disks of several exchange lengths in size. The energy of this large vortex state is computed numer- 
ically (within the framework of Magnetism@home distributed computing project) and its stability 
is studied analytically, which allows to plot it on magnetic phase diagram. It is the ground state 
of cylinders of certain sizes and is metastable in a wider set of geometries. Large vortices exist 
on par with classical ones, while being separated by an energy barrier, controllable by tuning the 
geometry and material of ferromagnetic disk. This state can be an excellent candidate for magnetic 
information storage not only because the resulting disk sizes are among the smallest, able to support 
magnetic vortices, but also because it is the closest to the classical vortex state of all other known 
metastable states of magnetic nano-cylinder, which implies, that the memory, based on switching 
between these two different types of magnetic vortices, may, potentially, achieve the highest possible 
rate of switching. 

PACS numbers: 75.60.Ch, 75.70.Kw, 85.70.Kh 
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Information storage is always based on switching some 
physical system between different metastable states, sep- 
arated by energy barrier. In the case of magnetic nano- 
cylinders with magnetic vortices (theoretically studied by 
N.A. Usov and S.E. Peschany[lj and confirmed experi- 



mentally by Shinjo et al[2 
vortex core polarizations 



the vortex states of different 
or different chiralitiesjj] are 
now being extensively researched for this purpose. There 
is much progress in this areaQ, including the remarkable 
results on switching the vortex core polarity by ultra- 
short in-plane magnetic field pulses [6 -8]. However, de- 
spite the driving pulses can be made extremely short, the 
switching itself is usually accompanied by much longer 
magnetization dynamics during which the new equilib- 
rium state is established If the switching process in- 
volves creation of additional vortices and anti- vortices, 
they must be annihilated at the end, which implies sig- 
nificant spin- wave generationfioj . The energy of these 
spin- waves (as well as other energy, accumulated by mag- 
netization) must be dissipated, which limits the sustained 
rate of switching of such a devices. This limitation is 
more pronounced, the longer is the trajectory, each spin 
must pass during the switching, or, the bigger is the dis- 
tance between the metastable states used. Knowing these 
states is, therefore, very important for applications and 
also is an intriguing fundamental problem. Similarity of 
the equations means that finding a new metastable state 
in nano-magnetism is like discovering a new elementary 
particle in non-linear field theory. With a difference that 
nano-magnets have boundary, allowing for existence of 
additional states, such as boundary-bound half- vortices 
and anti- vortices llj as well as large vortices. 



tures are non-linear and non-local, solving them for 
stationary states usually involves significant guesswork. 
Usually, new states (such as domain structures or new 
types of domain walls) are first observed in experiments 
(or numerical experiments) and only then described the- 
oretically. The experiment also does not provide re- 
searchers with a "silver bullet" , since there are many 
metastable states in magnetism and discovering a new 
one experimentally (often by accident) involves applying 
a certain (sometimes quite complex) sequence of external 
forces to the system (or, setting initial conditions in the 
case of numerical experiments) , which (may or may not) 
drive it to a new stable state. This task becomes espe- 
cially hard if the new state occupies only a small part of 
magnetic phase diagram or appears only as metastable 
(as opposed to the lowest energy ground state of the sys- 
tem). In case of nano-magnets, however, it is possible 
to parametrize the low-energy magnetization distribu- 
tions using several scalar parameters ll| , which facili- 
tates systematic exploration of their state space. Such 
exploration of states in circular cylinder was performed 
in the framework of distributed computing project Mag- 
netism@home. 

Specifically, the project is based on the following 
parametrization of magnetic texture: 
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Because the equations for equilibrium magnetic tex- 



which specifies components of dimensionless magneti- 
zation vector m(t) = M/Ms, normalized by satura- 
tion magnetization of cylinder's material Ms, via stere- 
ographic projection m x + im y = 2//(l + \f\ 2 ), m z — 
(1 — |/| 2 )/(1 + |/| 2 )- The texture is assumed to be uni- 



2 



form along the cylinder axis (Z) and depends only on 
the in-plane complex coordinate t = X + iY, where X, Y 
and Z are dimensionless Cartesian coordinates, normal- 
ized by the cylinder radius R so that \t\ — X 2 + Y 2 < 1. 
Along Z axis the cylinder is located at — g/2 < Z < g/2, 
where g = Lz / R is normalized thickness (aspect ratio) . 

The texture ([lj is aparticular single- vortex case of a 
more general texture H( with additional rescaling, gov- 
erned by the parameter p > 1 , to cover states with quasi- 
uniform magnetization[12j]. It can be applied to cylin- 
ders of arbitrary shape of their face, which is specified 
by the conformal map T(t) from unit disk |t| < 1 to 
the desired shape. In this paper (and in the first run of 
the Magnetism@home project) only the circular cylinders 
are considered, that is T(t) — t and T'(t) — 1. Complex 
phase of parameter a controls direction of the vortex cen- 
ter displacement. In circular cylinder this phase can be 
nullified by rotating the coordinate system. The parame- 
ter a was, therefore, assumed to be real < a < oo. The 
parameter < c < oo controls the size of the vortex core 
(or half- vortex cores, if they are present). 

Magnetism@home project numerically computed mag- 
netostatic energy (including that of volume, face and side 
magnetic charges, as well as their mutual interaction) of 
magnetization distributions (JT|) in a 4-dimensional unit- 
hypercube with normalized coordinates g — <?/(l + .<?), 
c = c/(l + c), a = a/(l + a), p = 1/p, which covers all the 
magnetization configurations it describes. The densities 
of magnetic charges were computed analytically, while 
the magnetostatic energy itself was evaluated using fast 
multipole method 13| on a dense 50000 non-uniform fi- 
nite elements mesh, perfectly covering the circular cylin- 
der's boundary, with exact analytical treatment for Z 
dependence of demagnetizing field. This computation 
(including preliminary testing runs) took about a year to 
complete, using idle processor time of a several tens of 
thousands of computers on the Internet, communicating 
via Berkeley Open Infrastructure for Network Comput- 
ing (BOINC) protocol. The exchange energy is much 
simpler to evaluate (and also its dependence on g can be 
taken into account analytically), which was done locally. 
The resulting data files (including the source code of pro- 
grams, used to compute and process them) are published 
as Magnetism® home data release 1, attached as supple- 
mental material. 

Combining and interpolating Magnetism@home files 
allows to compute multi-dimensional energy landscape 
of the particle as function of p, a and c in particles of 
different physical dimensions R/Le and Lz/ Le, where 
Le = y/C/^B^oMs is the exchange length of cylinder's 
material with exchange stiffness C, po is permeability of 
vacuum and 7b = 1 (in CGS system of units [l4| Ho = 1 
and 7_b = 47r). Then, building the ground state magnetic 
phase diagram is as simple as finding the smallest value in 
the resulting energy array at each point (R/ Le,Lz / Le) 
and classifying the corresponding state (given by values 



of p, a and c at minimum) . The resulting numerical dia- 
gram is shown in Figure by shading. 

Potentially, the files allow to extract much more in- 
formation, but for now, let us focus on the ground state 
diagram, as it already displays a new simple and fun- 
damental result: there is another magnetic vortex state, 
which, in tall cylinders (Lz > 6.5) of certain radii, is, 
actually, the ground state. At larger cylinder radii the 
energy of this state becomes equal (shown by the solid 
line a in Figure) and then larger than the energy of the 
classical magnetic vortex. At smaller radii it continu- 
ously transforms into the state of completely uniform 
out-of-plane magnetization (solid line b). This expands 
the region of existence of magnetic vortices to particles 
of smaller radii and solves a well known paradox 2 1(| 
that (when the large vortices are not taken into account) 
the uniform out-of-plane state continues to be the ground 
state of the cylinder (with energy, smaller than that of 
classical vortex) with geometries (below the dotted line 
c and above line 6), where it is no longer stable with re- 
spect to transformation into magnetic vortex. The line 
b corresponds to the second order phase transition dur- 
ing which the uniform state loses stability and total Z- 
component of magnetization starts to deviate from Ms- 
The line a corresponds to the first order transition, ac- 
companied by hysteresis, so that the large vortices may 
exist as metastable in certain region above it, while the 
classical vortices may exist below this line (they continue 
to be stable with respect to core expansion in cylinders 
with radii down to equilibrium classical vortex core ra- 
dius, shown in Figure as dash-dotted line py). 

To confirm the above numerical results and also to fully 
outline the region of stability of the large magnetic vor- 
tices let us now proceed with analytical consideration. 

Two interactions, present in any ferromagnet, are the 
exchange and magnetostatic self-interaction. In contin- 
uum approximation, assuming the magnetization texture 
is simply a function of coordinates m(r), r — {X,Y, Z}, 
the exchange energy can be represented as 
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(2) 



where integration is performed across the volume of the 
magnet V, outside of the magnet |m| = 0. In the consid- 
ered case of vortices, larger then the particle radius, the 
magnetization texture can be described by the following 
analytical function of complex variable 
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Pv 



(3) 



with py = Rv I R > 1 is the normalized vortex core ra- 
dius Ry- Since pv > 1, this magnetization texture is 
all-soliton[llj. That is, contrary to the case of classical 
magnetic vortex, there is no region in the particle, where 
magnetization vector lies in the cylinder plane. In this 
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FIG. 1. Equilibrium and stability of large magnetic vortex state in circular cylinders of different radii R and thicknesses Lz 
in soft magnetic material with exchange length Le and saturation magnetization Ms- Solid lines separate the regions (also 
having different texture), where different states (uniform in-plane, uniform out-of- plane, classical vortex and large vortex) have 
the lowest energy. They correspond to particle geometries, where the energy of respective pairs of states is equal. Cylinder 
geometries, where large magnetic vortex can be metastable, are shown by darker shading. Lettered lines, outlining this stability 
region, are discussed in the text, they correspond to different modes of vortex stability loss. The line p v shows the radius of 
classical vortex (which depends only on cylinder's thickness), computed by Usov and Peschanyp]]. Inset shows out-of-plane 
magnetization component distribution Mz/Ms as function of reduced radial coordinate r/R in large and classical vortex at a 
particular geometry, where they both are stable and have the same energy (E L = E c ). 



case, the exchange energy can be expressed directly in 
terms of f(z) as 
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which, after integration, gives 
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where the dimensionless energy e = E / (hq^bM^kLzR 2 ) 
and dimensionless particle radius p = R/Le were intro- 
duced. 

To compute the magnetostatic energy let us employ 
magnetic charges formalism, following from Maxwell's 
equations in the static case with no macroscopic currents 
present. Then, it is possible to express the demagnetizing 
field via gradient of a scalar potential, which is solution 
of Poisson's equations with the volume density of mag- 
netic charge f2 — — divM on the right hand side. On 
the boundary of magnetic material this volume density 
reduces to localized surface density of magnetic charge, 



equal to the component of magnetization vector, normal 
to the boundary. In the case of centered large vortex 
([3]) only the surface magnetic charge on the particle face, 
proportional to out-of-plane component of magnetization 
a = M z is present. Its energy can be expressed as 

cms = - (w{p v ,0) - w(p v ,g)), (6) 
9 



(5) where the magnetostatic function 
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(7) 

with r and ip being dimensionless coordinates on cylin- 
der's face. For centered vortex ^ the density of face 
charges does not depend on ip and the integral over the 
angles can be factored out, producing 
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where K(k) is complete elliptic integral of the first kind. 
The function w is a nice continuous function and can be 
differentiated over p v under the integral. 

To find equilibrium radius of large magnetic vortex it is 
now sufficient to differentiate the total energy e ex + e-MS 
over pv and require that this derivative is equal to 0, 
which leads to the following transcendental equation 

+ «— ( w (pv, o) - w [pv,g)) = o. (io) 



p 2 (i + p 2 v ) 2 dpv 

This equation is hard to solve for py. However, it is easy 
to solve for p, which gives 
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allowing to compute the radius of the cylinder p in which 
the large vortex of specified radius py would be at equi- 
librium. The limit py — > oo of this expression imme- 
diately recovers the line b in Figure - the boundary at 
which the large vortices (at larger radii) lose their chiral- 
ity and become the uniform out-of-plane state (at smaller 
radii) . 

The other stability boundary of the large vortex state, 
line d, can be found in two different ways. One is di- 
rect but cumbersome, by considering the transformation 
between the all-soliton magnetic vortex and the classical 
one, looking for the particle geometry where the former 
becomes unstable with respect to transformation into the 
latter. This can be achieved by considering magnetiza- 
tion distributions, expressed via the complex function 



w(t,t) = 
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(12) 



This function is, in general, not analytical. At a = 
it describes all-soliton vortex and at a = ir/2 non- 
analytical everywhere classical vortex, having a meron 
part at \t\ > py. The line d can be obtained as a bound- 
ary at which the energy minimum at a = disappears. 
Such analysis was indeed performed, however, it turns 
out (at least to numerical precision) that all-soliton vor- 
tices lose stability right at the moment their core bound- 
ary touches the particle boundary. Then they immedi- 
ately transform into classical vortices of much smaller 
core radius Py. Thus, the line d can be found simply by 
computing the particle radius (|lll) in the limit py — > 1. 

There is another catch that the line b we have just com- 
puted becomes unphysical when the aspect ratio of the 
particle g = Lz/R becomes smaller then approximately 
1.81295. At thinner cylinders the uniform in-plane state 
becomes the ground state and it becomes necessary to 
study the transformation of the large vortex state into 
it. This can be done by considering the uniform tilt of 
the magnetization in the large vortex, described by the 



complex function 
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icos(a/2) + /(f) sin(a/2)' 



(13) 



At a — it is a pure large vortex ([3]), but, as a increases, 
it's magnetization uniformly rotates towards the particle 
plane. The stability line, then, corresponds to the dis- 
appearance of the energy minimum at a — and can be 
found by solving the equation 
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at equilibrium py. Since the / tllt: (f) is still an analytical 
function of t the exchange energy e E \ can be directly 
obtained from ((4|). The magnetostatic energy, however, 
now has volume and side surface charge contribution in 
addition to the energy of face charges. Moreover, there 
is an interaction between the side and volume charges 
(their interaction with face charges is canceled by sym- 
metry), which also depend on the polar angle ip. This 
makes it necessary to adopt another method of comput- 
ing magnetostatic energy as angular integrals can't be 
directly factored-out. First, let us merge the volume and 
side charges density together, using Dirac's delta func- 
tion S(x) to express the side charges localization. Their 
product (entering the magnetostatic integral) is 
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The integral of such product of charges can be very effi- 
ciently evaluated (l7j by representing the inverse distance 
in polar coordinates via Lipshitz internal and then using 
the Bessel's summation theorem to factor the angular 
integral out, which gives 
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where J\(x) is Bessel's function of the first kind and the 
delta function is assumed to be right-sided. This integral 
was evaluated numerically, which is sufficient to solve the 
transcendental equation (fT4|) and compute the stability 
line e. Below this line the large vortices, if created, imme- 
diately convert into the uniform in-plane state. A confir- 
mation of the consistency of the performed computations 
is the fact that the line e joins the line b exactly at the 
critical aspect ratio g = 1.81295 at which the energies of 
uniform in-plane and out-of-plane states are equal. 

Finally, there is another possibility of vortex conver- 
sion into the uniform in-plane state, not captured by the 
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expression (|13p . At very small thicknesses the vortex 
may become unstable with respect to lateral shift. The 
stabilizing force in this case is produced by the side mag- 
netic charges, whose role quickly diminishes as particle 
becomes thinner. The corresponding stability line / can 
be found by equating to zero the second derivative in 
a at a = of the energy of magnetization distribution, 
corresponding to the following complex function 



r hift w 



Pv 



(18) 



The face charges energy is unchanged in this case, and 
the balance of forces is established between the exchange 
(|4| and magnetostatic interaction of side charges. The 
exchange prevails for radii below the stability line / and 
the particle becomes uniformly magnetized in-plane. 

Let us now discuss limitations of the above consider- 
ation. In reality there will be a deviation from the as- 
sumption of magnetization uniformity across cylinder's 
thickness. However, its influence on the phase dia- 
gram lines (in the considered range of thicknesses) is 
negligible, as evidenced by 3-d micromagnetic simula- 
tions and experiment (l8| . Also, the metastability re- 
gion of large vortex state comes down to Lz/Le ~ 0.3, 
where magnetization thickness independence is beyond 
any doubt. Full 3-d models can, however, bring in new 
effects. At large thicknesses the opposite faces of cylinder 
become magnetostatically decoupled, so in thicker cut 
cones (such as commonly produced by lithography), it 
might be possible to have a metastable large vortex state 
on smaller face, while having only the small classical-like 
vortex stable on larger face. Such states can be useful 
e.g. for readout of stored information, but their exis- 
tence needs a separate confirmation. Lithographic pro- 
cess may also produce other defects, such as imperfect 
boundary or local variations of thickness. While large 
deviations of this kind would produce completely unpre- 
dictable magnetic configurations, small random defects 
usually pin down states across the first order phase tran- 
sition boundary and make hysteresis loops wider. Ex- 
perimental points in Fig 14 of Ref. HH also indicate that 
dots with dimensions, supporting large vortices, are ex- 
perimentally accessible and are not superparamagnetic. 
The exact lifetime analysis of these states still needs to 
be performed, but from fundamental point of view they 
are viable and must be kinetically stable at least at suf- 
ficiently low temperatures. 

Concluding, it was shown that there is another vortex 
state in ferromagnetic circular cylinders with vortex core 
radius, larger than the particle radius. In certain range 
of particle geometries this large vortex may coexist with 
classical vortex [lj], whose core is fully inside the parti- 
cle. In this region one may have the energies of both 
classical vortex and large vortex equal (see inset in Fig- 



ure) while their magnetic moment differs significantly. It 
should be possible to switch magnetization between these 
two vortex states (e.g. by applying current through the 
particle). Moreover, because these states can be contin- 
uously transformed into one another without formation 
of Bloch points or other topological singularities and are 
otherwise very similar to each other, such switching can 
be expected to be very fast. 
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